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ABSTRACT 

We present chemical elemental abundances for 36, 561 stars observed by the RA- 
dial Velocity Experiment (RAVE), an ambitious spectroscopic survey of our Galaxy at 
Galactic latitudes |b| > 25° and with magnitudes in the range 9< -^denis <13. RAVE 
spectra cover the Ca-triplet region at 8410-8795A with resolving power R~7500. This 
first data release of the RAVE chemical catalogue is complementary to the third RAVE 
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data release of radial velocities and stellar parameters, and it contains chemical abun- 
dances for the elements Mg, Al, Si, Ca, Ti, Fe and Ni, with a mean error of ~0.2 dex, 
as judged from accuracy tests performed on synthetic and real spectra. Abundances are 
estimated through a dedicated processing pipeline in which the curve of growth of indi- 
vidual lines is obtained from a library of absorption-line equivalent widths to construct 
a model spectrum that is then matched to the observed spectrum via a x^-minimization 
technique. We plan to extend this pipeline to include estimates for other elements, such 
as oxygen and sulfur, in future data releases. 

Subject headings: catalogs — stars: abundances — techniques: spectroscopic — Galaxy: 
abundances — Galaxy: evolution — surveys 



Introduction 



Stars inherit the chemical patterns of the interstellar matter from which they were born. At 
the end of their lives they return their nuclear products to the interstellar medium through stellar 
winds and supernovae, enriching it with heavier elements. The elemental abundance pattern of 
every generation of stars thus depends on the previous one and, in time, a sort of "genealogical 
tree" develops. In principle, this allows for the star formati on his tory of a galaxy to be traced using 
stellar elemental abundances (Freeman &; Bland-Hawthorn |2002| ) . 



At the same time, the assembly history of a galaxy also leaves traces in the kinematics of its 
stars. Indeed, once born, stars behave like a collisioness fluid that spreads through phase space, 
generally leaving clear tracks of their dynamical origin. For example, disrupted open clusters 
generate moving groups of stars with similar kinematics that can still be recognized long after 
disruption (Eggen ll965l ). while accretion events pro duce t ransient stellar streams that closely track 
the orbit of their progenitor satellites (Helmi et al. , 120061 ) . 



Chemical and kinematic information, when available, can therefore be used to reconstruct the 
history of the Milky Way much in the way archaeologists examine relics to recreate the ancient 
past. The data requirements of this exercise, however, have limited its applicability in the past: as- 
sembling a large and homogeneous set of proper motions, distances, radial velocities, and elemental 
abundances requires large photometric and spectroscopic surveys that have not been feasible until 
recently. 

Spectroscopy has traditionally been the bottleneck, with much of the work on chemical elemen- 
tal abundances limited to small, biased samples. Until recently, iron abundances ([Fe/H]) had been 
estim ated f or samples as large as the ~ 16, 000 stars of the Geneva-Copenhagen Survey (Nordstrom 
et al.,|200J), but the numbers of stars with abundances measured for multiple elements were much 
smaller. Inde ed, the largest homogeneous sample available up to now was published by Valenti & 
Fisher (j2005l ). who measured the abundances of 5 elements in 1,040 nearby F, G and K stars. 
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Smaller samples, often compos ed of a few hund red sta.rs or fewer, have also been prese nted by 
several auth ors (Edvardsson et al.. ll993l . Puhrmann [l99i, ISQO^, Luck & Heiter bood . bopj Reddy 



et al. |2006| . among others). La rger, but inhomogeneou s, dat asets have been collated from the 



literature by Soubiran &: Girard (j2005l ) and Venn et al. (|200J) , who compiled 743 and 821 stars, 
respectively, and also in the PASTEL catalogue (Soubiran et al.. l201ol ). which is a collection of 865 
literature studies. 

The availability of dedicated telescopes and multi-object spectrographs has radically changed 
thi s stat e of affairs, enabling surveys such as the RAdial Velocity Experiment, RAVE (Steinmetz at 
al. l2006l). an d the Sloan Extension for Galactic Understanding and Exploration (SEGUE) (Yanny 
et al., |2009|). Combined these two surveys have now taken spectra for roughly a million stars. 
Galactic archaeology has thus become one of the fastest-growing fields of astronomical enquiry, 
as evidenced by the nume rous surveys currently underway or in the advanced planning stages, 
e.g., Gaia (Ferryman at al. l200ll ): the Large Sky Area Mul ti-Obj ect Fiber Spectroscopic Telescope 



(LAMOST0); and HERMES (Freeman & Bland-Hawthorn I2H) . 



We present here the first release of the chemical catalogue for RAVE, an ambitious spec- 
troscopic survey of our Galaxy at Galactic latitudes |b| > 25° and magnitudes in the range 
9< -^j^ENjg <13. This catalogue contains multi-element abundance measurements for 36,561 
stars of the Milky Way based on 37, 848 RAVE spectra covering the Ca-triplet region at 8410- 
8795A with resolving powe r R^ 7500. This data release is associated to the RAVE third data 
release (DR3, Siebert et al. I2OIII ). where further information (kinematics, photometry, etc.) can 
be found. As of summer 2011, RAVE has taken more than half a million spectra for some 400,000 
stars. Abundances for these targets will be published in subsequent data releases. 



RAVE (Steinmetz et al. |2006| ) was first conceived as a radial velocity survey to provide the 
missing third velocity component for stars in the solar suburb. However, it soon became clear 
that RAVE spectra carry much more information than just radial velocities. After further develop- 
ment, the RAVE processing pipeline was modified to deliver estimations of the value s of the stellar 
param eters like effective temperature, gravity and metallicity (Zwitter et al. l2008l . Siebert et al. 



2011 



)• 



We take this development one step further here by tackling the measurement of chemical 
elemental abundances at the medium spectral resolution provided by RAVE data. This feature 
is shared with other large spectroscopic surveys with limited resolution and spectral coverage, 
such as Gaia and LAMOST. Although high precision is not expected from medium resolution 
spectroscopy, the availability of hundreds of thousands of RAVE spectra enables the creation of a 
large and homogeneous catalogue of chemical abundances suitable for statistical investigation. 

The RAVE chemical pipeline measures abundances for 7 elements: Mg, Al, Si, Ca, Ti, Fe 
and Ni. Homogeneity is assured by an automated processing pipeline that measures abundances 



""^ http://www.lamost.org/website/en 
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assuming the stellar parameters computed by the RAVE pipeline described in Sie bert et al. (|201ll ). 
RAVE stars have astrometry fro m diff erent sources like Tycho2 (H0g et al. l200d ). PPM- Extended 
catalogue PPMX (Roeser et al. 200^ and the Second U.S. Naval Observatory CCD Astrograph 
Catalog UCAC2 (Zachari as et al. 120041 ). The se, to gether with R AVE r adial velocities and distance 
estimates (Breddels et al. 2010 . Zwitter et al. 2010 . Burnett et al. boill ). yield 3D positions and ve- 
locities. Combining this information with chemical abundances results in a unique chemo-kinematic 
dataset suitable for investigations of the formation history of the Galaxy. 

This paper is structured as follows. In Section [2] we present the chemical processing pipeline, 
detailing how the chemical elemental abundances are measured. Tests to establish the accuracy 
and reliability of the results are detailed in Section [3l while in Section [5] we discuss our method of 
measurement and in Section [5] we present the RAVE chemical catalog. In Section [6] we outline our 
conclusions. 



2. The Pipeline 

The RAVE chemical abundance pipeline uses a different approach to classical elemental abun- 
dance estimation methods. These methods either (i) measure equivalent widths (EW)s and infer 
elemental abundances from the curves-of-growth (COG)s or (ii) synthesize spectra with varying 
elemental abundances to find the best match between the synthetic and observed spectra. The 
first method cannot be successfully applied to the RAVE spectra because it requires isolated lines 
for precise EW measurements and, at the medium resolution of RAVE, most of the lines are in- 
strumentally blended. Also, the second method is computationally too expensive as it requires the 
synthesis of several spectra ~400A wide, each having hundreds of lines. 

The adopted method can be considered a hybrid approach, measuring the elemental abun- 
dances by fitting the spectrum with a model constructed with lines of known EWs. The construction 
of the model is robust and fast, with the elemental abundances determined by a minimization 
routine. The drawback of this method is that the EWs of the lines are computed by neglecting 
the opacity of the neighboring lines. This may overestimate the EWs, leading in some cases to 
underestimated elemental abundances. We show below that this effect is small for most of the lines 
in the RAVE wavelength range. 



2.1. Procedure 

The chemical pipeline utilizes RAVE reduced spectra and their stellar parameters, such as 
effective temperature (T^^^^), gravity (logg'^^^^) and metallicity ([m/H]^^^^), returned by the 
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RAVE pipeline (described in Siebert et al., I2OIII 1 I. The pipeline also uses an EW library based 



on the RAVE line list. The RAVE line list contains 604 known absorption lines identified in the 
RAVE wavelength range (Sec l2.2p . The EW library contains the EWs of absorption lines for every 
point on the T^g, log 5 and [m/H] grid (Sec. 12. 3p . which accounts for the variation of the EW with 
the stellar parameters. Variations of EW with elemental abundances are then further calculated 
for each point of this grid, for five different abundance levels in the range [—0.4, 0.4] with respect 
to the value of metallicity [m/H]. 

The pipeline contains some auxiliary codes to (re) normalize the spectrum and detect spectral 
defects like bad normalization, cosmic rays and other unwanted features. 

The pipeline algorithm can be outlined as follows: 

1. Upload the normalized, RV corrected and wavelength calibrated spectrum and the estimated 
stellar parameters T^'J^'"', log^"^^'^, [m/H]"^^'^ 

2. Upload the RAVE line list from the EW library and the corresponding EWs for the stellar 
parameters at the five different abundance levels 

3. Extract a shorter line list of those lines which, at the estimated stellar parameters, have large 
enough EWs to be visible above the noise 

4. Fit the strong lines and correct the continuum 

5. Construct the COGs of the lines by fitting a polynomial function through the five EW- 
abundance points 

6. Create the model by assuming a Gaussian profile for each line and summing these profiles 
together 

7. Vary the chemical elemental abundances to obtain different models by changing the EWs of 
the lines according to their COG 

8. Minimize the between the models and the observed spectrum to find the best-matching 
model 

In the following we give further details pertaining to these steps. 



2.2. Line list and constraining the oscillator strength 

The RAVE line list contains 604 absorption lines identified in spectra of the Sun and Arcturus. 
The lines correspond to the element species N I, O I, Mg I, Al I, Si I, S I, Ca I, Ti I, Ti II, Cr I, Fe I, 



[m/H]'''^^^indicates the uncalibrated metallicity as defined in Zwitter et al. 20081 . 
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Fe II, Co I, Ni I, Zr I and to the CN molecule. In order to get precise chemical elemental abundances 
from the EWs we firstly need reliable atomic parameters for the lines. A critical parameter is the 
oscillator strength (gf, often expressed as logarithm log gf). If precise laboratory measurements 
are missing for the log gf values then they are obtained through an inverse spectral analysis to 
obtain 'astrophysical lo g q/s'. This is the case for the RAVE wavelength region, where the VALD 



database (Kupka et al. Il999l ) reports that only 11 Fe Ilines have laboratory oscillator strength 
measurements. The Sun and Arcturus's spectra were synthesized and the log gf values for each 
line obtained using an automated procedure which varies the log gfs so as to match both spectra 
simult aneously. For the synthesis we adopted Solar abundances as given by Grevesse & Sauval 



(jl998l ). They will be the zero point of the abundances of this work. The details of the procedure 



are outlined in Appendix lAl 

The pipeline makes use of a shorter working line list which contains the lines strong enough 
to be visible above the noise. The lines deemed visible are those satisfying: 

EW{mA)>2.^.d-1000 

where STN is the signal-to-noise raticjfl and d is the spectral dispersion (for RAVE spectra d = 
0.4A/pix). Other lines can blend with these visible lines, contributing to their flux absorption and 
changing their apparent EW. We therefore added to the working line list all those lines having 
EW>lmA and being closer than 0.6A (1.2A is the full-width half maximum of a typical RAVE 
spectrum) to the visible lines. 



2.3. The EW library and COG computation 

In order to compute the line COGs we built an EW library with the RAVE line list for 
30, 640 atmos phere models. The EWs were computed usi ng th e LTE line-analysis software MOOG 



(Sneden ll973l ) and the ATLAS9 model atmosphere grid (|2003l . Castelli & Kurucz). As the spacing 
of the ATLAS9 grid is different to our requirements, we linearly interpolated the ATLAS9 grid to 
produce our own grid, covering the T^a range [3600,7600]K in steps of lOOK, gravity range [0.0,5.0] 
dex in steps of 0.2 dex, and metallicities in the range [—2.5,-1-0.5] dex in steps of 0.1 dex. We did 
not compute the EWs for atmospheres with Tcfl>5100K and log(7<1.0 (hot supergiants), because 
MOOG does not produce results for these parameters and furthermore there are no such stars in 
our catalogue. For every atmosphere model we computed the EWs of the lines for 5 abundance 
levels with respect to iron: [X/Fe]= -0.4, -0.2, 0.0, 0.2, 0.4 dex assuming [Fe/H] = [m/H]. The 
whole EW library consists of 145, 080 files. To obtain EWs for stellar parameters between the grid 
points of the EW library we linearly interpolate the closest points on the grid. 



''in the following we will use STN to indicate the signal-to-noise ratio as computed in Siebert et al. (201l|) whereas 
S/N will be used with the usual meaning given in spectrophotometry. 
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Fig. 1. — Distribution on {T^g, log g) plane of 712 stars observed spectroscopically at high resolution. 



2.4. Microturbulence 



The EW computation requires a value of the microturbulence for each atmosphere model. 
For high-resolution data is typically determined by measuring the EWs of Fe Ilines and changing 
until the iron abundances inferred from strong and weak lines agree. As we do not measure 
individual EWs for our spectra we cannot use this procedure, but instead rely on relations giving 
£ as a function of th e stell ar parameters. Such re lation s have been derived by Edvardsson et al. 
(jl993l ). Reddy et al., (j2003l ). Allende Prieto et al., (j2004l ). where ^ is given as a function of T^g and 
log 5. Unfortunately, these results cover only specific regions of parameter space (e.g., hot dwarfs or 
cold giants). Thus, we derived our own relation covering a wide a range of T^g and log g. To do so we 
made use of literature results from high resolution spectroscopy that report both stellar parameters 
and £ of their stellar samples. We col lected data for 712 giants and dw arfs spanni ng a wide range 
in Teff and log g from Luck k He iter bood . l2007l . Bensby et al. boosl . Fuhrmann llQod . Fulbright 
et al. l2006l . Allende Prieto et al. 12004 . Figure [1] displays the coverage in the T^a - log 5 plane for 
the sample. A 3rd-degree polynomial fit was used to obtain the microturbulence dependence on 
gravity and effective temperature. Appendix [B] gives the coefficients of this polynomial fit, (,poiy 



Figure [3 compares the Allende Prieto et al. (j2004l ) law with (,poiy Uncertainties of cr^=0.32 
km s~^ in our polynomial law translate approximately i nto ~ 0.04 dex elemental abundance un- 



certainties for dwarfs stars (as estimated by Reddy et al. 120031 and Mishenina et al. 120031 ) 



2.5. Continuum re-normalization and strong line fitting 



Before measuring EWs, the spectrum's continuum needs to be re-normalized and its strong 
lines fitted. For while the spectrum is continuum normalized by the RAVE pipeline, we have found 
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that this normahzation is is not rigorous enough for elemental abundance estimation: the RAVE 
pipeline employs a 3rd order spline function to fit the continuum leaving behind fringing effects 
on scales of less than 50A wide seen in some RAVE-normalized spectra. The chemical pipeline 
performs a new normalization to remove such unwanted features, fitting the continuum and the 
strong Ca Hand H iPaschen lines. These strong lines are then excluded from the measurement 
process because they are difficult to fit properly. The pipeline then considers as "continuum" the 
fit to the classical continuum plus the strong lines, and by comparison with this level the metallic 
lines are measured. 

Before applying the continuum correction we estimate a preliminary metallicity, which we 
call [m/H]°''' to distinguish it from the final metallicity [m/H]'''"'™. This preliminary metallicity is 
required because the intensity of the absorption lines must be known in order to subtract them. 
When they are subtracted the continuum level can then be seen. However, in order to measure 
the intensity of a line one must know where the continuum lies. The continuum level is therefore 
determined iteratively, starting with the normalization of the RAVE pipeline and measuring the 
line intensities using the chemical elemental abundances estimation subroutine (see Sec. 12. 6p . with 
the difference that all the abundances vary together as one variable [X/H]=[m/H]°°*. Wavelength 
intervals ~20A wide centered on the strong lines are avoided as the large wings of these lines can 
affect the results of the determination of the metallicity. The best-match [m/H]'''" is then used 
to synthesize the metallic lines, which are then subtracted before the strong lines are fit. 

Fitting the continuum is performed in four steps, summarized below and illustrated in FigureO 
from top to bottom: 

1. Perform a preliminary metallicity and subtract the metallic lines (gray line in panel (a). 
Figure Ej) . 

2. Fit the strong lines with a Lorentzian profile for Ca Hand a Gaussian profile for H land 
subtract them (gray line in panel (b). Figure [3|). 

3. Estimate the continuum profile by box-car smoothing what remains (gray line in panel (c). 
Figure ED. 

4. Add the strong lines and the continuum together to obtain the new "continuum" (grey line 
in panel (d). Figure [3]) 



2.6. Chemical elemental abundance estimation 

The chemical pipeline uses the parameters T^^^ and logg^^^^ for the elemental abundance 
estimation. The observed spectrum is fit by a model spectrum firstly obtained by subtracting 
the fiux absorbed by the metallic lines to that derived above in Sec. 12.51 the strong lines 
-|- continuum fit (gray line in panel (d). Figure [3]). The absorption lines are assumed to have 
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Gaussian profiles, and as the instrumental profile is dominant with respect to the line profile, we 
use the same full width half maximum (FWHM) for all lines. However, this FWHM is varied and 
optimized during the minimization process. 

The pipeline computes the COG of the lines by using the EWs at five levels of abundances 
([X/Fe]= -0.4, -0.2, 0.0, 0.2, 0.4 dex, assuming [Fe/H] = [m/H]^^^'=) from the EW library. For 
every line, the five EW points are fit with a 3rd order polynomial function, which serves as the 
COG. However, these polynomial functions represent only part of the COGs and will diverge if 
extrapolated too far beyond the five levels in the EW library. To avoid such divergences we limit 
accepted abundance results to those in the range —0.6 <[X/Fe]< +0.fE|. 

Using the COGs the pipeline then creates the model spectrum, using [X/H] = [m/H]'^^^^ as a 
first-guess to the metallicity. It computes the between the observed spectrum and the model 
and, through a minimization process, changes the elemental abundances [X/H] until the best match 
(minimum x^) with the observed spectrum is reached. The minimization process is performed with 
15 variables: 13 elemental abundances, one molecule and the instrumental FWHM of the lines. 

Figure |4] shows the best-match model spectrum of the Sun and compares it with the observed 
one. The three spectra on the top show how the model spectrum is built up: the spectra of three 
elements (Fe, Si, Mg) at parameters T^g =5861K and logg' =4.54 are generated according to the 
estimated elemental abundances and added together to construct the solar spectrum. 

This methodology has a drawback: the computed COG of any absoption line neglects the 
opacity of the neighboring lines, leading to an overestimation of the EW of blended lines and so 
an underestimation of the elemental abundances. This systematic originates from the EW library 
created by using the MOOG's driver "ewfind", which computes the expected EW of the lines as if 
they were isolated. In order to minimize this systematic error we apply a correction coefficient to 
reduce the EWs of lines which are physically blended. The corrected EW is given by 

EWcorr = EW ■ coeff ■ cont, (1) 

where the coefficient 

neighbor<0.2A 

coeff = 1 — EWi/2.50/dispersion, 

i 

is computed by considering all the neighbouring lines within 0.2A of the line of interest. The 
dispersion is expressed in A/pix and cont is the value of the continuum (gray line in panel (d) of 
Figure [3]) at the central wavelength of the line. The multiplication of the EW by cont in Eq. [1] 
corrects for the effects of strong lines (such as Ca II), if the line of interest lies within their large 
wings. This correction reduces the systematic error in the abundances from ~ —0.15 dex to ~ —0.1 
dex or less for most elements (see details in the following sections, where quality checks are outlined). 



^Abundances beyond this limit can be due to photon noise stronger than the absorption hues or by pecuhar 
abundances. The two cases cannot be distinguished by the chemical pipeline. 
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2.7. Spectrum quality determination: the frac parameter 

Roughly 25% of RAVE spectra are affected by defects like fringing or cosmic rays that cannot 
be removed by continuum correction (they usually affect regions smaller than lOOA). In order to 
determine the locus and size of the defect we define the flux residual between the observed and the 
best-match model for the i-th pixel with 

r{i) = fmodelii) - fobsii), 

wher e / is the flux of the spectra, and we use the following algorithm (also used in Siebert et al., 
201lh : 

1. Consider the interval Ij = [j 

^ max{Ij) 

V'(j) = 77T V r(i), 

^ ^' ^ i=rain{lj) 

where x^(i) is the reduced and ^(j) is the estimation of the area between the observed 
and model spectra, a is the inverse of the signal-to-noise (STN) ratio. 

2. If x^(j) > 2 and V'(j) > 2 • a then the pixels in the interval Ij are labeled as a defect. 

The process is repeated for all the pixels of the spectrum, resulting in an array of and 1, 
with 1 label indicating non-defective pixels. The fraction of pixels that are non-defective is then 
denoted with the jrac parameter (i.e., the higher the frac parameter, the better the quality of the 
spectrum). We deem spectra with frac < 0.7 as overly affected by fringing/cosmic ray defects and 
exclude these spectra from our analysis. 



10, j + 10] centered on the pixel j. Compute 

max{Ij) 

\ Al) ^ 

max(Ii) — min(Ii) ^-^ a ' 

^ ■'^ ^ ■'^ i=min(Ij) 



3. Validation and accuracy 



There are several sources of error that affect the elemental abundances measurement process. 
For example, there are uncertainties in (i) the oscillator strength values, (ii) the continuum nor- 
malization, (iii) the stellar parameters; (iv) uncertainties in microturbulence, (v) uncertainties due 
to local line opacity neglected by the pipeline, (vi) neglect of non-LTE effects and (vii) photon 
shot noise. These sources of errors interact in a complex fashion, making the error estimation pro- 
cess quite challenging. To establish the capability of the pipeline to derive the chemical elemental 
abundances we ran a series of tests on synthetic and real spectra for which stellar parameters and 
elemental abundances are well known. The results are presented below. 



3.1. Chemical elemental abundances accuracy from synthetic spectra 

We ran the pipeline on a sample of 1,353 synthetic spectra to which three intensity levels of 
artificial noise were added, testing the accuracy of the results at S/N= 100, 40, and 20. To make 
the sample as realistic as possible, the spectra have been synthesized with distributions of T^^ and 
log (7 from a mock sample of RAVE obs ervati ons (M. Williams, private communication) created by 



using the Besangon model (Robin et al. l2003l ). The T^^ vs log 5 distribution of the sample is shown 
in Figure El 

Each spectrum of the sa mple had the elemental abundances of one of the stars from the 



Soubiran & Girard catalogue (j2005l ). ensuring realistic distribution of the chemical abundances. 
As our line list has only eight elements in common with the Soubiran &: Girard catalogue, we 
assigned the elemental abundance [X/H] = [Fe/H] to the elements Cr, Co, Ni, Zr. We also assigned 
[X/H] = [a/H] to C, N, O, Mg, Al, Si, S, Ca, Ti when the measurement of one of these elements 
was missing. The spectra were synthesized using the code MOOG at resolution O.OlA/pix and 
degraded to RAVE resolution (0.4A/pix, 1.4A FWHM). 

The first test was performed by adopting the stellar parameters T^^, log 5 and [m/H] of the 
synthetic spectrum in order to evaluate the accuracy of the elemental abundance measurements of 
the pipeline alone. The results are illustrated in Figure[6l For the second test we added uncertainties 
to the stellar parameters so as to simulate real RAVE spectra. Following the typical RAVE errors, 
we added a random Gaussian error with standard deviation 300K in T^ff, 0.5 dex in log (7 and 0.3 
dex in [m/H] to the correct parameters and ran the pipeline adopting such "wrong" parameters. 
For the sake of brevity and clarity we present only the results for the seven most reliable elements 
(see also Figure [7| . 



3.1.1. Results assuming no error in the stellar parameters 
Results at S/N=100. With no errors added to the stellar parameters the residuals between 
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measured and expected elemental abundances [X/H] have a systematic error of ~ —0.10 dex and 
standard deviation cr ~ 0.10 dex or smaller (see Figure [6|), with a slight variation from element to 
element. The relative enhancement [X/Fe] shows a small (~ —0.10 dex) systematic error at high 
[X/Fe]. 

Results at S/N=40. The elemental abundance residuals have a ~0. 10-0.20 dex depending 
on the element and show the same systematics seen with S/N=100, but slightly less pronounced. 

Results at S/N=20. At this S/N the pipeline can still estimate abundances of Fe, Al and 
Si and the relative enhancement to iron [X/Fe], but with a ~0.2-0.3 dex. Other elements like Mg 
and Ti exhibit larger systematic errors. An a-element abundance estimated as the average of Mg 
and Si still yields reliable results with an error of o" ~ 0.2 dex. Elements such as Ca and Ni cannot 
be reliably measured because the lines are too weak to be detected. 

We found correlations between stellar parameters and elemental chemical abundances obtained. 
In particular, the abundances correlate with T^ff (Figure [8]). This is very likely due to continuum 
correction, which can appear lower than the real level in spectra crowded of lines, as in cold giants 
stars. 

3.1.2. Results including errors in the stellar parameters 

Results at S/N=100. When errors are added to the stellar parameters the elemental abun- 
dances [X/H] are affected to the level of a ~ 0.15 — 0.20 dex for elements like Mg, Al, Si, Ca, Fe, 
Ni. Titanium abundances show clear systematic errors, being overestimated at high abundance 
and underestimated at low abundance. Enhancements relative to Fe, [X/Fe], have errors a < 0.2 
dex and little sign of systematic errors for most of the elements. 

Results at S/N=40. The elemental abundance errors are in the range a ~0. 15-0.30 dex 
depending on the element. The relative elemental abundances [X/Fe] are reliable with errors 
cj ~0.2 dex. 

Results at S/N=20. At this S/N the pipeline can still estimate abundances of Fe, Al and 
Si and the relative enhancement with respect to iron [X/Fe] even with an error of a ~0.2-0.3 dex. 
Other elements like Mg and Ti show systematic errors. An a-element abundance computed by 
averaging Mg and Si yields reliable results with an error of o" ~ 0.2 dex. Elements such as Ca and 
Ni cannot be reliably measured because the lines are too weak to be detected. 

As could be expected, there is a correlation between errors in stellar parameters and errors in 
the elemental abundance estimates. Figure [9] highlights these correlations, showing that most of 
the elemental abundances errors are correlated with T^g errors. In general, an overestimate in T^^ 
corresponds to an overestimate in abundance, since most of the lines have smaller EW at higher 
T^ff. Only weak correlations with log 5 and [m/H] errors are evident. 
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3.2. The accuracy of elemental abundance estimates from real spectra 

The initial RAVE input catalogue had no stars whose chemical elemental abundances were 
known to high accuracy to compare the results of our pipeline with. To create a co mpari son dataset 



of real stars, we therefore observed 104 stars chosen from Soubiran & Girard (120051 ) (hereafter 
SG05), a collection of high-precision elemental abundance measurements from the literature. Since 
these are bright stars, most of the RAVE spectra have S/N>100 and therefore the a ccuracy in 



our abundance estimates is representative of the high S/N case. Ruchti et al. (|2010l ) (hereafter 
RIO) have also re-observed 243 RAVE stars at high resolution and measured stellar parameters 
and chemical elemental abundances. The RAVE spectra of these stars have S/N between 30 and 
90 and therefore we can use them to test the accuracy of our procedure for intermediate S/N. 

We have seven elements in common with SG05 and five elements with RIO (Fe, Mg, Si, Ca and 
Ti). However, for the latter we could only compare four elements because the weak Ca Ilines were 
not strong enough in the RAVE spectra to be measured by the pipeline. As with synthetic spectra, 
we ran a first test by adopting the high-resolution stellar parameters to evaluate the performance 
of the pipeline alone and then ran a second test by using the RAVE stellar parameters, which have 
larger uncertainties relative to the high-resolution data. 

The results of the first test are shown in Figure [10] where the gray dots and black "+" corre- 
spond to SG05 and RIO, respectively. When the high-resolution stellar parameters are adopted the 
chemical pipeline's elemental abundances agree with the SG05 and RIO results to within ~ 0.15 
dex on average. The smaller dispersion of the SG05 stars around the one-to-one correspondence 
line compared to RIO is due to the higher S/N of the RAVE spectra of these stars. All elements are 
typically underestimated by ~0.1 dex, an offset that is roughly constant across the entire abundance 
range [—2.0,-1-0.5] dex. The a-enhancement (Figure [TT]) is estimated well for all but [a/Fe]>0.4 
dex, where there is a small systematic bias of ~0.1 dex. 

When the RAVE stellar parameters are adopted (the "second test" , Figure [T2]) the dispersion 
in our elemental abundances compared to the published values increases to ~0.2-0.3 dex on average, 
depending on the element. In Table [U we report the mean and standard deviation after grouping 
the sample in dwarf and giant stars separately. The larger dispersion at low abundances ([X/H]<- 
1.0 dex) is very likely due to the fact that estimating precise stellar parameters on such spectra 
is challenging due to the few and weak visible absorption features (excluding for the Ca Iltriplet, 
which is not measured). 

The enhancement estimates ([X/Fe], see Figure [13]) are accurate to within ~0.3 dex. At high 
enhancements, abundances are systematically underestimated by ~0. 1-0.2 dex, depending on the 
element. However, as we are comparing our measurements with other authors' (high resolution) 
measurements, the variance of the residuals is the quadratic sum of the variances of our and the high 
resolution results. This means that the standard deviations reported in Figure [12] and Figure [13] 
represent a conservative estimate of the expected RAVE errors. 
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[Si/H] - 
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0.21 
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0.14 
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0.30 


[Fe/H] - 
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0.15 
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0.30 


0.04 


0.31 


[Ni/H] - 
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-0.02 


0.29 






0.06 


0.26 


[a/H] - 


[a/H]ex 


-0.10 


0.15 


-0.04 


0.10 


-0.02 


0.21 


-0.02 


0.15 



Table 1: Mean and standard deviation of the residuals between measured and expected elemental 
abundances obtained by using the high resolution stellar parameters (left) and the RAVE param- 
eters (right) for 347 standard stars. We computed the values for giants (log(7<3.5) and dwarfs 
(logg'>3.5) separately. The two groups are also representative of stars with T;,ff <5500K and inter- 
mediate STN (RIO stars), and stars with T,«>5500K and high STN (SG05 stars). 



3.3. Internal errors from repeat observations 

In order to take into account the effect of binary stars, RAVE observes some stars (~5% of 
the whole sample) more than once. In this chemical catalogue, 964 stars have been observed twice, 
90 thrice and 67 four times. Multiple observations are also useful to estimate the internal error 
by comparing the results obtained from different spectra of the same object. A fair comparison 
would require that the repeated observations of an object have the same S/N. Unfortunately, due to 
variable atmospheric conditions the spectra often have different S/N values. We therefore averaged 
the STNs of the spectra for each object and computed the corresponding standard deviations of 
the elemental abundances. In this way we could assess the magnitude of internal errors in our 
procedure as a function of STN. In Figure [l5] we plot the results for 12, 504 RAVE stars having 
multiple observations, where we use the RAVE internal data release as it has a substantially larger 
number of stars. For all the elements the la confidence line (gray solid line) lies below 0.2 dex and 
hovers around ~0.1 dex for STN>80. 



3.4. [m/H]^'^^'= vs. [m/H]""^""": a comparison 



In Figure [m we compare the distributions of [m/H]'^^^'', [m/H]"'^""" and [Fe/H]'^'"™^ for 37819 
RAVE spectra with STN>20. The metallicity [m/H] i s inferred from the chemical elemental 
abundances with the equation given by Salaris et al. (119931 ) 

[Fe/H] + log(0.638 • lol"/^^' 



[m/H] 



chem 



0.362), 



(2) 
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where the a enhancement is computed as 

KFel^»Wa_pe/H, ,3, 

and the elemental abundances [Mg/H], [Si/H] and [Fe/H] come from the chemical pipeline. 

[m/H]^^^^ is ~0.1 dex lower than [m/H]''''™% which is to be expected as the non-calibrated 
RAV E meta llicity is underestimated by ~0.15 dex with respect to the reference stars used in Siebert 



et al.. l201ll . The shape of the [m/H]^*^^ distribution is fairly similar to the [m/H]''''°™ distribution 
for dwarf stars but different for giants. Moreover, [m/H]^*^^ seems to better match [Fe/H]"''™ than 
[m/H]'''"'™, particularly for giants. This could be due to a enhancement: giants stars have a higher 
proportion of thick disk stars (a-enhanced) than do dwarfs stars, which are mostly thin disk (not 
a-enhanced). Although used during the stellar parameters estimation process, the RAVE pipeline 



is unable to measure a enhancement (Zwitter et al. . 120081 ). Therefore, it is possible that particularly 
for a-enhanced stars [m/H]^^^^ better represents [Fe/H]. Apart from the discussed shift of 0.1 dex 
in metallicity, there seems to be fair agreement between the [m/H]^^^^ and [m/H]'''"'™ distributions. 



4. Discussion 

For each element, the accuracy of the abundance estimate depends on the number of lines that 
are strong enough to be clearly detected; to be measurable, the lines must be stronger than the 
noise. At the same time, the intensity of the lines depends on the stellar parameters and elemental 
abundances of the star. This means that for a fixed S/N, the accuracy of the elemental abundance 
is a function of the stellar parameters and of the abundance itself. 

The interplay of these factors makes the accuracy of the abundance results difficult to estimate. 
At high S/N, the accuracy varies from element to element, depending on the number and intensity 
of the lines measured. For instance, Ti and Ni lines are measured better in stars with T^g lower 
than -^SOOOK, whereas sulfur (not present in this data release) can be seen and measured only 
for temperatures higher than ~4800K. The RAVE line list has a few weak Ca Ilines that can be 
measured only in stars with [m/H] >-0.5 dex. 

Noise affects the EWs of the lines and so the model spectrum. As the pipeline varies the 
abundances of the individual elements, the larger the number the absoption lines of an element, 
the lower the abundance error due to the fitting, because the effects of noise are averaged over the 
lines. Thus, the abundance accuracy of elements with fewer lines is especially affected by noise. 
When the noise is particularly strong with respect to the intensity of the lines, i.e., low S/N, the 
minimization routine fits the noise and the resulting elemental abundance diverges. In such cases 
the pipeline rejects the measurement and renders a null value -9.99. 

The situation get worse in the low S/N regime, because the pipeline measures even fewer 
lines, i.e., only those strong enough to be detected despite the noise. This decreases the number of 
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measurable elements and increases the uncertainties. In Section [3. II we tested the pipeline down to 
S/N=20 to see if at such a low S/N the measurements are still trustworthy. The results suggest we 
can use such data, but with care. Noise generates selection effects: spectra having low [m/H] or 
high T^ff do not have lines strong enough to overcome the noise and they go through the pipeline 
unmeasured. A limited number of lines can also generate systematic errors because the abundance 
of an element depends on just a few lines (sometimes even only one) and thus becomes especially 
vulnerable to uncertainties in the oscillator strengths. 

Nevertheless, the tests performed with synthetic spectra at S/N=20 show that abundances of 
elements with strong lines like Fe, Al and Si can still be measured. At such low S/N, errors as large 
as ~0.2-0.3 dex are expected, but they do not show large systematic biases and the residuals are 
on average close to zero. 

We conclude that elemental abundances may be trusted down to S/N=40 for seven elements 
(details will follow in Section [5]), whereas between S/N=20 and 40 we can trust the abundance of 
[Fe/H] and (to be on the safe side) the abundance [a/H], when computed from the average of Mg 
and Si. Abundances for other elements for the low S/N stars should be considered indicative, and 
used only for exploratory purposes until new comparison stars enable a proper validation of our 
results in this S/N regime. 

4.1. Zero point of the RAVE elemental abundance scale 

Our elemental abundance measurements are indirect measurements in the sense that they 
are inferred from the comparison between the intensity of lines seen in real spectra and their 
intensity predicted by stellar atmosphere models. Since the models are different for different stellar 
parameters, a question arises regarding whether all the models yield elemental abundances which 
refer to the same zero point, i.e., the origin of the internal elemental abundance scale. The same 
question concerns whether this zero point refers to the same zero point of the real spectra, i.e., 
between the internal and external scales. 

The latter question has a prompt answer: we do not know the external scale, because real 
stellar atmospheres have never been directly probed and all elemental abundance measurements 
refer to models. Therefore, we can only check the consistency of the internal scale. This can be 
performed by comparing the measured elemental abundances of a sample of synthetic spectra at 
different stellar parameters, as done in Figure [6l 

In this plot, at S/N=100, the residuals between measured and expected elemental abundances 
are on average zero for any gravity and metallicity; the points align along a straight line with slope 
roughly equal to one. This means that the measured differences in elemental abundance at any 
metallicity regime are the same (constant offset), i.e., they refer to the same zero point. However, 
the offset from element to element may vary, showing that the offsets are due to the measurement 
process and not due to the metallicity of the atmosphere models. 
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On the other hand, there is a clear correlation between residuals and T^s. the higher the T^ft ,the 
higher the measured elemental abundance. This systematic may be due to continuum correction: 
cooler stars have spectra that are crowded with absorption lines and so the continuum appears 
lower than it is. Additionally, the varying behaviour of the elements in Figure [6] can be explained 
by considering that their lines lie in different regions of the spectrum, i.e., on the wings of a Ca Illine 
or in a region relatively free of lines. This can translate into systematic, wavelength-localized shifts 
of the estimated continuum as a function of T^g. 

Although this systematic effect appear tractable and correctable, we did not apply any correc- 
tion to the data because such a systematic is not visible on tests with real spectra. When the RAVE 
elemental abundances are compared with the expected elemental abundances of standard stars as 
function of T^ff no clear trend is visible (Figure [T4l) , except that RAVE elemental abundances have 
a slight systematic bias to lower values. A reasonable explanation is that the elemental abundance 
measurements given by SG05 and RIO suffer from the same systematic error, due to uncertainties 
in estimating the correct continuum level in regions crowded with absorption lines or affected by 
the wings of strong, broad lines. 

We conclude that the RAVE chemical pipeline can determine chemical elemental abundances 
with a systematic error that may be as large as ~ 0.1 dex as a function of T^g. Since this systematic 
looks quite linear with T^^ (see Figure [9] and Figure [T4|l . to reduce this error to ±0.05 dex we 
suggest analysing separately samples of stars with T<,fi>5500K (mostly dwarfs) from samples with 
T^ff <5500K (mostly giants). 

5. The RAVE chemical catalogue 

We present the catalogue of chemical elemental abundances for 36,561 RAVE stars, measured 
using 37,848 spectra. 

5.1. Sample selection 

The spectra have been selected from the DR3 RAVE database using the following constraints: 

• Effective temperature 4000<Te,f(K)<7000: this is the temperature range within which 
the RAVE line list has been calibrated. At lower temperature the spectra are characterized by 
molecular lines other than CN (CH, TiO and other molecules), while at higher temperatures 
lines of ionized atoms appear, and both are not included in the line list. 

• Signal-to-noise STN>20: for STN<20 the absorption lines are strongly affected by noise 
and the stellar parameters and chemical elemental abundances are not reliable. 
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• Rotational velocity V^ot <50 km s : at higher rotational velocity the lines have a FWHM 
larger than that due to the spectral resolution (FWHM~1.2A) and they cannot be precisely 
measured. Moreover, any spectra showing larger lines might be a double-lined spectrum, 
which we also must be avoided. 

Despite these selection criteria, some spectra exhibit emission or are affected by bad continuum 
normalization. For such spectra the stellar parameters and chemical elemental abundances are not 
reliable and it is advisable to reject them. Constraints on the parameters ^-^d frac help identify 
these spectra. For statistical studies we suggest rejecting spectra with > 2000 and frac < 0.7, 
which removes 739 out of 37,848 spectra of the catalogue. 



5.2. Stellar parameters 

The distributions of the stellar parameters and STN values are presented in Figure [T71 The 
quantities T^^ and log (7 are given by the RAVE data archive whereas STN and [Fe/H] are estimated 
by the chemical pipeline. The distributions in T^g and log g show two peaks corresponding to giants 
at lower temperature and dwarfs, which are mostly at higher temperatures. The iron abundance 
distribution peaks at [Fe/H]~-0.1 dex for dwarfs and [Fe/H]~-0.5 dex for giants. The latter are 
in average more metal-poor because they lie further from the Galactic plane and have therefore a 
larger fraction of thick disk starjfj. 



5.3. Chemical elemental abundances: selection effects due to STN and stellar 

parameters 

We encourage the user to pay particular attention to selection effects introduced by the stellar 
parameters and STN. Together they affect the total number of abundance estimations as well as 
their accuracy. Absorption lines can be undetectable because of the low [m/H] of the star, because 
the too-high (or too-low) T^^ does not populate enough of its electronic levels, or because the 
spectrum has a too-low STN. This is illustrated in Fig [18] and Fig 1191 where the number of spectra 
having [X/H] estimation diminishes with STN and [m/H]. In general, only metal-rich stars have 
elemental abundance estimations at any STN whereas metal poor stars have estimations only if 
their spectra have high STN. 



^As RAVE is a magnitude limited survey, the intrinsically bright objects are on average at large distances. 
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5.4. Accuracy and reliability element by element 

We now discuss and summarize the reliability of the chemical abundances for individual ele- 
ments in light of the previous discussion. First, we make some general remarks about the measured 
elemental abundances and their errors: 

• The accuracy of the chemical abundances depends on several variables. In particular if T^g is 
erroneously estimated, the abundances are affected to different degree for different elements. 
If there are systematics deviations in T^^ there will be systematic deviations in [X/H] as well. 

• In general [X/H] are underestimated, and the magnitude of the bias is a function of T^h' 
stars with T^g <5000K yield on average abundances that are '^0.1 dex lower than stars with 
>5000K, whose abundance errors average to zero. This effect increases the errors when stars 
with different T^g are analyzed. On the other hand, relative abundances [X/Fe] are nearly 
unaffected because the trend is similar for all elements. 

• The errors given below refer to the expected errors for two intervals of STN at [m/H] ~0.0 
dex. As already discussed in the previous section, errors can increase for lower STN and lower 
[m/H]. 

Magnesium yields reliable results on synthetic and real spectra. At any T^g we expect an abun- 
dance error aug <0.15 dex for STN>40 and ~0.25 dex for 20>STN>40. 

Aluminum abundances are obtained from only two physically isolated lines (which are instrumen- 
tally blended at RAVE resolution). Despite the instrumental blending, the lines are strong and 
give accurate estimates in tests with synthetic and real spectra that shows no systematic offset. 
For STN>40 we expect abundance errors ctai <0.2 dex and ~0.3 dex for 20>STN>40. 

Silicon is the most reliable element together with Fe. We expect an abundance error crsi <0.15 
dex for STN>40 and ~0.20 dex for 20>STN>40. 

Calcium abundances arc obtained from only five weak Ca Ilines. Ca lis better measured at higher 
metallicity and T^g <5000K. Estimated errors are (Xca ~0.25 dex for STN>40 and ~0.4 dex for 
20>STN>40. 

Titanium gives reliable abundances at high STN. At STN~40 its abundance is reliable for Teff<5000K 
and underestimated for higher temperatures. The correlation with T^g errors is particularly strong 
(Teff underestimation generates [Ti/H] underestimation and vice versa), leading to larger errors. 
We expect an abundance error of crxi ~0.2 dex at STN>40 and ~0.3 dex for 20>STN>40. We 
recommend separately analyzing stars with T^g lower and higher than 5000K. 

Iron is the most reliable element together with Si. It can be accurately measured on spectra with 
any T^g in the range wc consider. We expect an abundance error of (jpe ~0.1 dex at STN>40 and 
-0.2 dex at 20>STN>40. 
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Nickel has six weak lines in the RAVE wavelength range that are visible only at r;,ff <5000K. In 
synthetic spectra with STN=100 the pipeline yields errors comparable to those of other elements 
but biased to lower values by ~0.2 dex. It is not measurable for [m/H]<-0.6 dex. Tests on real 
spectra are inconclusive because they were all performed for stars with T;,fi>5000K. We expect 
uncertainties of dNi ~0.2 dex for STN>40 and ~0.3 dex for 20>STN>40. Because of the smah 
number and the weakness of its lines, it is advisable to use Ni values with care. 

The overall metal licity [m/H]''''''™ is inferred from the chemical elemental abundances with Equa- 



tion[2](Salaris et al. ll993l ). with the abundances [Mg/H], [Si/H] and [Fe/H] delivered by the chemical 
pipeline. The metallicity [m/H]'''""" is reliable and we expect errors of Ufj^/H] ~0.1 dex for STN>40 
and ^0.2 dex for 20>STN>40. 

The a enhancement [a/Fe] has been computed as the average of [Mg/Fe] and [Si/Fe] (Equa- 
tion [3]). In the range 20>STN>40 it is advisable to use it instead of the single a element values 
because it is more reliable. We estimate errors of a a '^0.1 dex for STN>40 and ~0.2 dex for 
20>STN>40. 



5.5. The data 

The RAVE chemical catalogue contains data for 37,848 spectra of 36,561 stars and it is provided 
as an ASCII table of 37,848 lines. It contains chemical abundances for the elements Mg, Al, Si, Ca, 
Ti, Fe and Ni, RAVE stellar parameters, signal-to-noise STN, object name and other parameters 
for quality checks as explained in Table [2j The distributions of the relative chemical abundances 
are shown in Figure [20l The catalogue will be accessible online at |http : //www . rave- survey . orgj 
and via the Strasbourg Astronomical Data Center (CDS) service. 
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Fig. 2. — Top panel: residuals between measured microturbulence ^ and computed ^ with the 
AUende Prieto et al.'s formula for 712 stars. Bottom panel: residuals between measured ^ and 
computed by a 3rd degree polynomial function. Solid and open points represent dwarf and giant 
stars, respectively. The gray cross represents the Sun. 
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Fig. 3. — Black line: RAVE spectrum of the Sun. Panel (a): metallic lines fitting. Panel (b): 
strong line fitting. Panel (c): continuum fitting. Panel (d): the grey line is the sum of the gray 
lines in panels (b) and (c). This curve is used as the continuum line from which we estimate the 
chemical elemental abundances. 
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Fig. 4. — Once the abundances of Fe, Si and Mg are fixed their absorption fines can be synthesized 

(gray fines on the top). By summing them together we obtain the grey fine on the bottom, which 
is the model spectrum of the Sun. The black line is the Sun's spectrum reflected by the Moon and 
observed by RAVE. 
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Fig. 5. — Distribution of the synthetic spectra sample on the stellar parameters plane T^g and log 5. 
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Fig. 6. — Left: expected elemental abundances [X/H] (x-axis) versus measured elemental abun- 
dances (y-axis) for the sample of synthetic spectra at S/N=100, 40, 20 (for the left, middle and right 
column, respectively) and assuming no errors in stellar parameters. Right: as in left panels but 
for the expected enhancement [X/Fe] (x-axis) and the residuals measured-minus-expected (y-axis). 
Offsets and standard deviations are reported in the panels. 
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Fig. 7, 
of the 



. — As in Figure [6] but with noisy stehar parameters to simulate the RAVE stellar parameters 
RAVE archive. 
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Fig. 8. — Correlation between the elemental abundance residuals measured-minus-expected (y-axis) 
and the stellar parameters (x-axis) at S/N=100. 
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Fig. 9. — Correlation between the elemental abundance residuals measured-minus-expected (y-axis) 
and the stellar parameter errors (x-axis) at S/N=100. 
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Fig. 10. — Comparison between high resolution elemental abundance ([X/H]j,g£) from the hterature 
(x-axis) and RAVE elemental abundances (y-axis) for the SG05 sample (98 stars, gray dots) and 
the RIO sample (243 stars, black "+"). For these measurements we adopted the values of T^a and 
logfjr from the high-resolution data (first test). 
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Fig. 11. — Relative elemental abundance residuals RAVE-minus-high resolution (y-axis) versus 
high resolution measurements ([X/Fe]ref ). For these measurements we adopted the values of T^g 
and log 5 from the high-resolution data (first test). The symbols are as in Figure [TOl 
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Fig. 12. — As in Figure [10] but adopting T^ff and \ogg from the RAVE data (second test). 
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Fig. 13. — As in Figure [TT] but adopting T^ff and logg from the RAVE data (second test). 
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Fig. 14. — Residuals between measured and expected elemental abundances as a function of T^a. 
For these measurements we adopted high resolution stellar parameters. 
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Fig. 15. — Expected internal errors of the chemical elemental abundances, estimated via standard 
deviations of the abundances from repeat observations. For a more robust statistic we used 30821 
spectra of 12504 stars of the RAVE internal data release. The solid gray line represents the la 
abundance error limit, i.e., 68% of the stars lie under the line for each STN bin, whereas the dashed 
gray line represents the 2a limit. For this sample 9469 stars have two observations, 1286 have three, 
1060 have four, 473 have five, 153 have six, 33 have seven, 11 have eight, 11 have nine, 4 have ten 
and four have 11 observations. 
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Fig. 16.— Comparison between [m/H]"'^""" (black line) [Fe/H]"'^""" (dotted line) and [m/H]^^^^ (non- 
calibrated metallicity, gray line) distributions for the whole sample (37848 spectra, left panel), dwarf 
(middle panel) and giant stars (right panel). 




Fig. 17. — Distributions of stellar parameters T^g, logg, [Fe/H] and signal-to-noise ratio STN for 
37848 spectra of the chemical catalogue. In the bottom left panel, the [Fe/H] distribution is given 
separately for dwarfs (black line) and giants (gray line). 
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Fig. 18. — Fraction of spectra having elemental abundance measurements as a function of the STN. 
Every bin is normalized to the total number of observed spectra in the corresponding STN bin. 




Fig. 19. — Fraction of spectra having elemental abundance measurements as a function of 
[m/H]^^^^. Every bin is normalized to the total number of observed spectra in that [m/H]"^^^^ 
bin. The solid, dashed and dotted lines represent the whole sample, dwarf stars only and giants 
stars only, respectively. 
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Fig. 20. — Abundances relative to iron for spectra with STN>20 (left panels) and STN>40 (right 
panels). The isocontours hold 34.0%, 68.0%, 95.0% and 99.5% of the sample. 
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6. Conclusions 

This first release of the RAVE chemical catalogue reports the chemical abundances of seven 
elements (Mg, Al, Si, Ca, Ti, Fe, Ni) for 36,561 stars (from 37,848 spectra) selected from the DR3 
data release. The chemical processing pipeline developed for the creation of the present catalogue 
assures homogeneous elemental abundance measurements in addition to the RAVE pipeline that 
estimates the values of the stellar parameters and radial velocities. 

The uncertainties in the estimated elemental abundances depend mostly on the signal-to-noise 
ratio (STN) of the spectra; on the values of the stellar parameters of the stars; and on the element 
considered. The expected errors in abundance can be as small as 0.1 dex at STN>80 for the 
elements Fe and Si and up to ~0.3 at 20>STN>40 for Ni. 

On average we expect errors of ~0.2 dex at STN>40 for most of the elements. The relative 
abundances [X/Fe] are reliable with an underestimation of ~ —0.1 dex for [X/Fe]>+0.3. 

We plan to increase the number of measured elements in future data releases. We believe that 
with further development of the chemical pipeline it should be possible to enrich the catalogue with 
the abundances of at least two other elements (oxygen and sulfur). 

The availability of chemical elemental abundances, together with the radial velocities of the 



st ars (gi ven in the RAVEDR3 data release) and their distances (Breddels et al. l20ld . Zwitter et 



al. I2OIOI . Burnett et al. l201ll ). make RAVE a prime dataset for Galactic archaeology investigations. 
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Table 2: Catalogue description 



Field number 


name 


null value 


description 


i 


vJDject iJJ 




Identifier of the star 




z 


[Mg/HJ 


Q QQ 

-y.yy 


Mg abundance 


q 


1\T 




number of Mg lines measured 


4 


[Ai/MJ 


Q QQ 

-y.yy 


Al abundance 


r 



1\T 
IN 




number of Al lines measured 


/? 



rc; /ul 

[Si/HJ 


-9.99 


Si abundance 


7 


TV 
IN 




number of Si lines measured 


Q 

o 


[Ua/rlJ 


n no 

-y.yy 


Ca abundance 


9 


N 




number of Ca lines measured 


10 


[Ti/H] 


-9.99 


Ti abundance 


11 


N 




number of Ti lines measured 


12 


[Fe/H] 


-9.99 


Fe abundance 


13 


N 




number of Fe lines measured 


14 


[Ni/H] 


-9.99 


Ni abundance 


15 


N 




number of Ni lines measured 


16 






RAVE effective temperature 


17 


log 5 




RAVE gravity 


18 






RAVE metallicity 


19 


[m/H] 




metallicity from the chemical pipeline 


20 


[ct/Fe]'^'"''"' 




a-enhanccment from the chemical pipeline 


21 


STN 




signal-to-noisc ratio 


22 


frac 




fraction of spectrum matching the template well 


23 


Ntot 




total number of lines measured 


24 






between observed and template spectra 



- 35 - 



A. Correction of the atomic oscillator strengths in the RAVE/Gaia wavelength 

region 

In the following we outline how the pipeline treats oscillator strengths {gf, often expressed as 
a logarithm loggf) in the wavelength range 8400-8800A. These are necessary in order to obtain 
trustworthy elemental abundance measurements in a wav elengt h region which is still poorly known. 
According to the Vienna Atomic Line Database VALD (119991 ) only 1 1 Fe lines have accur ate q f - 



value laboratory measurements (Blackwell et al., 119861 . Bard et al., Il99ll . Bard & Kock, 11994 ). 
while the remaining hundreds of lines rely on theoretical and semi-empirical calculations (Kurucz, 
19951 ). Some of these values have proved imp recise , with errors as large as 1 dex for important lines 
belonging to Fe and Si (Bigot & Thevenin, |2006| ). We have attempted to improve the oscillator 
strength values of 604 absorption lines using an inverse spectral analysis applied to spectra of the 
Sun and Arcturus. 



A.l. Method 



The RAVE line list consists of 604 absorption lines (excluding the H land Ca I lstrong lines) 
identified by looking at the high S/N Sun and Arcturus s pectra (Hinkle et al., l2000l) and verified 
by cross checking the line identification of Moore et al. (jl966l ) and Hinkle et al. (j200d ). Some 
previously unidentified lines have been added in order to match features that are clearly visible in 
spectra of the Sun and Arcturus. 



The lines have been selected from Kurucz data (Kurucz, Il995l ) and we h ave also adopted 
Kurucz's wavelengths and excitation potentials x- Following Chen et al. (j200ol ) we also adopted 
the enhancement factor E^ for the Unsold approximation because it is commonly accepted that 
the line broadening due to Van der Waals interaction obtained from Unsold's approximation is too 
weak and needs to be corrected. For the CN molecule we used the dissociation energy D=7.75eV, 
as suggested by Sauval & Grevesse (119941 ). 



We performed the inverse spectral analysis by synthesizing sol ar and Arcturus's spectra with 
MOOG (a standard LTE line analysis and synthesis code by Sneden. ll973l ) and changing the loggf s 
values until a good match with both observed spectra is reached. We show here that matching these 
two spectra (or more) allows us to improve the loggf s of blended lines, which otherwise would be 
impossible to disentangle. 

The EW of a measur ed abs orption line is related to the loggf value through the curve of 
growth equation (see Gray, |2005| . formula 16.4) 



log 



EW 



log C + log ^ + log{gf ■ X)-ex- log Ku, 



(Al) 
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where C and are functionqj of the temperature T, and is the opacity at the frequency v. When 
all the atmosphere model parameters are fixed, the curve of growth equation can be written as a 
function, /, of log gf and elemental abundance A 

EW = /{log gf, A). (A2) 



This allows us to disentangle two or more blended lines: let EW^^^^^ be the EW of a blended fea- 
ture of the star x composed of n lines of m different elements (re > m) and having EWi, EW2---EWn 
unknown EWs. Then we can write the system of equations 

n 

EW^lend = j;/(log5/n,^n) (A3) 

n 

EW^lend = 5;/(log5/«,^n) 
n 

EW^lend = 5]/(log5/«,^n) 



where A^ represents the abundance of the element the line belongs to. 

Since the physical conditions and elemental abundance A of each stars are in general different, 
the values EW^i^^^ are different as well. The equation system [A3] is solvable when the (A'') number 
of equations is greater than the unknown re + rre. The existence and uniqueness of the solution is 
assured even for lines where the weak-line approximation does not hold, because the first derivative 
of the curve of growth is always positive. 

This means that it is always possible to determine the loggfs of a blended feature and the 
elemental abundances by using a large enough number of spectra of different stars. In our case we 
have used two spectra: the Sun and Arcturus. 



A. 2. The correction routine 

The following routine is valid when the stellar parameters and elemental abundances of both 
stars are known. We also consider later the case where elemental abundances are not known. To 



C and 9 are expressed as follows 

u{T) 



C = constant ■ ^^4^N„ 



9 = 5040/r 

where T is the temperature, Nj/Ns is the fraction of atoms of the j-th stage of ionization with respect to the number 
of atoms of the element considered, Nh is the number of hydrogen atoms per unit volume and M(r) is the partition 
function. 
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synthesize the solar spectrum we assumed an effective temperature Toff=5777K, gravity log ^=4. 44, 
microturbulence Vf=1.0 km/s and metalli city [m /H]=0.00. For Arcturus the atmospheric parame- 



ters are taken from Earle, Luck & Heiter (|2005|) (r,ff=4340K, logc/=1.93, Vt=1.87 km/s [m/H]=- 
0.55 dex). 

Th e synt hetic spectra have been obtained assuming solar chemical abundances by Grevesse & 
Sauval |l998|) (hereafter GS98) and adopting the grid ATLAS9 model atmospheres of Castelli & 
Kurucz (|2003i ). The loggfs correction is performed using the following procedure: 

1. Synthesize the solar and Arcturus spectra from our line list. 

2. Select the first line of the list 

3. Integrate the absorbed flux of the observed and synthesized spectra over an interval spanning 
O.4A around the center of the line (roughly 1 FWHM, which takes into account most of the 
absorbed flux) and calculate the normalized EW residual (NEWR), defined as follows 

NEWR = ~ ^"'^^ (A4) 

-^obs 

where Fgynt and Fobs ai'e the flux of the synthetic and the observed spectrum, respectively. 
For the Sun and Arcturus we write (NEWRsun) and (NEWRArc)- 

4. If NEWRsun < -0.05 then add +0.01 to the log 5/, else if NEWRsun > 0.05 then add 
—0.01 to the log gf (i.e. match the Sun first). Go to step 7. 

5. If \NEWRsun\ < 0.05 and NEWRatc < -0.05 then add -0.01 to the log 5/, else if 
\NEWRsun\ < 0.05 and NEWRatc > 0.05 then add +0.01 to the log gf (e.g. if the Sun 
matches then try to match Arcturus). Go to step 7. 

6. If \NEWRsun\ and \NEWRa rc\ are both < 0.05 then a valid log gf has been reached. 

7. Select the next line from the line list and return to step 3. 



The routine is iterated until no more improvement is seen. For about ten lines a good match 
could not be reached. Several reasons may be behind this: bad continuum correction, lines located 
over the Ca Illines where the fit is unaccurate, misidentified lines or unrecognized blends. For 
these lines, we did a visual check and manually corrected the log gf values. The correction is good 
enough to derive gf-values that yield abundance errors per measured line smaller than 0.2 dex in 
average for both the Sun and Arcturus. 

Since the abundance of Arcturus is not known for all elements, we decided to assume an initial 
abundance for all the elements ([X/H] = [m/H]) and find the right abundances with the following 
procedure: 



1. Synthesize the Sun's and Arcturus's spectra. 
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2. Correct the log gf values by running the above correction routine until no more improvements 
are possible. 

3. Check the distribution of the NEWRs for each element separately: if a positive/negative 
offset is present, decrease/increase its abundance. 

4. Synthesize the solar and Arcturus's spectra with the new loggfs and abundances. 

5. Return to step 2. 

This algorithm is repeated until there are no more improvements. The results of this procedure are 
the new loggfs and Arcturus's chemical abundances (listed in Table [3]). 



A.2.1. Multiplets treatment 

In the RAVE/Gaia wavelength range there are features that belong to atomic multiplets. An 
interesting characteristic of these multiplets is that they share the same excitation potential x and 
have small differences in wavelength. In this case, obtaining a proper log gf for the individual lines 
is problematic. 

We work around this problem in the following way: consider the two Al Ilines at Ai = 
8773.8961 and A2 = 8773. 897A that share the same excitation potential x = 4.021947el/. The two 
lines have unknown equivalent widths EWi and EW2 but the sum 

EWTot = EWi + EW2, (A5) 
is measured and well known. We can write the curve of growth as 

EW C-Aai- {gf) ■ A 



A 



(A6) 



where A^i is the aluminum abundance and the other quantities are as in Eq. lAll Substituting 
Eq. IA6I in Eq. IA5I and adopting the approximation Ai ~ A2 = A we obtain 

EWTot _ EWi + EW2 _C- Am ■ {gfi + 5/2) • A 



A A lO^x . 

Using the logarithmic form, this equation then becomes 



(A7) 



EWTnt 

log = log C + log A + log((5/i + 5/2) • A) - - log K., (A8) 

which suggests that an appropriate functional form for a "dummy" variable loggfs would be 
log(g/i + 5/2)- Even if this does not allow us to disentangle the two loggfs, we can consider 
the multiplet as if it were one single line and correct the log gfjic, it will be a "dummy" value but 
it results in giving the correct abundance of Al. This correction has also been applied to three 
Mg Imultiplets, centered at wavelengths 8717.814, 8736.020 and 8773.896A. 
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abundance 


this work 


LH05 


T98 


F07 


[0/H] 


+0.17 




-0.30 


-0.10 


[Mg/H] 


-0.03 


+0.02 


-0.45 


-0.10 


[Al/H] 


-0.27 


-0.28 


-0.40 


-0.17 


[Si/H] 


-0.15 


-0.14 


-0.30 


-0.15 


[S/H] 


-0.13 








[Ca/H] 


-0.26 


-0.56 


-0.20 


-0.28 


[Ti/H] 


-0.22 


-0.39 


-0.20 


-0.20 


[Cr/H] 


-0.37 


-0.55 


-0.20 




[Fe/H] 


-0.52 


-0.55 


-0.40 


-0.54 


[Co/H] 


-0.17 


-0.36 


-0.20 




[Ni/H] 


-0.45 


-0.48 


-0.35 




[Zr/H] 


-1.00 









Table 3: Arcturus' s elem ental abundances comp ared with the results by Luck & Heiter (j2005l ) 
(LH05), Thevenin (|l998l ^ (T98), Fulbright et al. (|2007l l (F07) 



A. 3. Tests 



In order to test the accuracy of this procedure, we have measured by hand the EWs of 100 
atomic lines in the solar spectrum. Figure [21] presents the difference in abundances obtained by 
using the Kurucz ^/-values (open points) and the corrected ones (solid points). Our correction 
reduced the standard deviations of the abundance residuals from cr=0.63 dex to cr=0.15 dex and 
their averages from +0.16 to -0.05. 



We also compared our results to the log gf by Kirby et al. (j2008l ). who corrected the ^/-values 
by applying a similar method, i.e., by synthesizing the solar and Arcturus's spectra and comparing 
the synthetic spectra with the observed ones. A comparison between our and their results is shown 
in Figure [22] for the 89 corrected lines in common: on average our log gf are lower by -0.05 dex. 
This is very likely due to different solar elemental abundances adopted: Kirby adop ts the Anders 
& Grevesse (Il989l ) abundances which are on average higher than Grevesse & Sauval (Il998l ). which 
we adopt. 

One more test involves comparing our results for Arcturus's elemental abundances with abun- 
dances quoted in the literature. This is do ne in Tab. [3] wh ere our elemental abundanc es are 
compared with those of Luck & Heiter (I2OO5I 1. Thevenin (Il998l ^ (T98), and Fulbright et al. ^20oi ) 
(F07). From the consistency of these tests, we are confident that our computed log gf values are 
accurate and that our procedure yields robust and reliable elemental abundance estimates. 
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EW Sun (mA) 

Fig. 21. — solar elemental abundances obtained from 100 atomic lines for which EWs (x-axis) have 
been measured by hand from the solar spectrum. The open points represent the elemental abun- 
dances obtained by using the Kurucz ^/-values, whereas the corrected ^/-values are represented 
by the solid points. The respective dispersions are (7=0.63 dex (open points) and (7=0.15 dex (solid 
points) . 




-4 -2 
log gfKirby 



Fig. 22. — Comparison between corrected loggfs by Kirby et al. 2008 (x-axis) and RAVE (y-axis). 
The statistics on the panel refers to points after the crossed outlier has been removed. 
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Table 4. Seven lines extracted from the RAVE line list. The whole table contains 604 lines with 

log gf correction, and the 3 Ca Hand 8 H Ilines which are not used for elemental abundance 
measurements (and did not undergo to log gf corrections). The atomic species 607.0 represent the 

CN molecule, following the MOOG format. 



wavelength(A) 


atomic species 


x{eV) 


log 9/ 




Dissociation energy (eV) 


8643.306 


26.0 


4.913304 


-2.470 


1.40 


0.00 


8644.077 


607.0 


0.933782 


-1.773 


2.50 


7.75 


8646.050 


607.0 


1.324726 


-1.540 


2.50 


7.75 


8646.371 


14.0 


6.206423 


-1.740 


1.30 


0.00 


8647.175 


607.0 


0.877831 


-1.872 


2.50 


7.75 


8648.465 


14.0 


6.206423 


0.100 


1.30 


0.00 


8648.688 


607.0 


0.910915 


-1.516 


2.50 


7.75 



Note. — Table |4] is published in its entirety in the electronic edition of the Astronomical 
Journal. 
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B. Microturbulence 

The polynomial function used to determine the microturbulence ^ is expressed as follows 

i+j<3 

^polyi^m s-^) = Yl ^ijT\loggy (Bl) 

ij=0,l,2,3 

where the coefficients 

aoo = -f 0.3533 (B2) 

aoi = 2.59492 

ao2 = 0.161863 

ao3 = 0.176579 

aio = 0.00509193 

ail = -0.00157151 

ai2 = -3.99782 • 10"^ 

= -0.000361822 

021 = 3.82802 • lO"'^ 

ago = -5.18845 ■ 10"^^ 
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